Draft version January 4, 2011 

Preprint typeset using L£-T^X style cmulateapj v. 11/10/09 



THE REDSHIFT SEARCH RECEIVER 3 MM WAVELENGTH SPECTRA OF 10 GALAXIES 

Ronald L. Snell 1 , Gopal Narayanan 1 , Min S. Yun 1 , Mark Heyer 1 , Aeree Chung 1,2 , William M. Irvine 1,3 , Neal R. 

Erickson 1 , and Guilin Liu 1 

Draft version January 4, 2011 

ABSTRACT 

The 3 mm wavelength spectra of 10 galaxies have been obtained at the Five College Radio As- 
tronomy Observatory using a new, very broadband receiver and spectrometer, called the Redshift 
Search Receiver (RSR). The RSR has an instantaneous bandwidth of 37 GHz covering frequencies 
from 74 to 111 GHz, and has a spectral resolution of 31 MHz (~100 km s _1 ). During tests of the 
RSR on the FCRAO 14 m telescope the complete 3 mm spectra of the central regions of NGC 253, 
Maffei 2, NGC1068, IC 342, M82, NGC 3079, NGC 3690, NGC 4258, Arp 220 and NGC 6240 were 
obtained. Within the wavelength band covered by the RSR, 20 spectral lines from 14 different atomic 
and molecular species were detected. Based on simultaneous fits to the spectrum of each galaxy, a 
number of key molecular line ratios are derived. A simple model which assumes the emission arises 
from an ensemble of Milky Way-like Giant Nolecular Cloud cores can adequately fit the observed 
line ratios using molecular abundances based on Galactic molecular cloud cores. Variations seen in 
some line ratios, such as 13 CO/HCN and HCO+/HCN, can be explained if the mean density of the 
molecular gas varies from galaxy to galaxy. However, NGC 3690, NGC 4258 and NGC 6240 show 
very large HCO + /HCN ratios and require significant abundance enhancement of HCO + over HCN, 
possible due to the proximity to active galactic nucleus activity. Finally, the mass of dense molecular 
gas is estimated and we infer that 25-85% of the total molecular gas in the central regions of these 
galaxies must have densities greater than 10 4 cm" 3 . 

Subject headings: galaxies: ISM — ISM: Molecules — ISM: Abundances — Radio Lines: galaxies 



1. INTRODUCTION 

The molecular emission lines observed in galaxies are 
powerful probes of the physical and chemical properties 
of the gas most directly connected to the star forma- 
tion process. Over the past several decades there have 
been many studies of the molecular line emission from 
galaxies, and these studies have provided invaluable in- 
formation on the molecular gas conte nt, star formation 
efficiency and mole cular abundances ([Young fc Scovillel 
119911: lOmontl 120071 ). Although CO emission has been 
the primary means of deriving the molecular gas content 
in galaxies, emission from other less abundant molecular 
species, is much better in determining the physical and 
chemical properties of the gas. As the chemistry of the 
molecular gas is affected by the local radiation field, these 
trace molecular constituents may also provide informa- 
tion on local environments of this gas. More frequently 
observations of molecular species such as HCN , HCO + , 
HNC, and CS are being obtained (lOmont| | 2007T) and in a 
few c ases mapped in nearby galaxies (IMeier fc Turner! 
2005). Since these molecules have much larger per- 
manent electric dipole moments than CO, they require 
considerable higher densities to collisionally excite, and 
thus, in general, trace denser molecular gas than the gas 
probed by CO. Emissions from these dense gas tracers, 
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such as HCN, have been argued to be a much better 
probe of the star-forming mole cular gas than emission 
from CO (|Gao & Solomonl[200l . 

The utility of emission from these 'high dipole moment 
molecules' to probe the properties of the nuclear regions 
of galaxies is well established, and there have been nu- 
merous papers on this subject in recent years (see re- 
view bv lOmontl ([20 07)). In addition, several emission 
line ratios, such as those for the low-lying transitions 
of HCN/CO, HCO+/HCN, and HNC/HCN, have been 
suggested as good diagnostics of t he properti es of the 
dense star-forming gas in galaxies ([Aaltol 120081 ). Nearly 
all studies of the molecular emission in galaxies have been 
limited to jus t a few targeted molecular lines. However, 
iMartin et alJ (|2006l ) recently presented the first spectral 
scan of a galaxy providing an inventory of the molecu- 
lar lines in NGC 253 within the 2 mm wavelength band. 
Such spectral scans can provide a much more complete 
description of the chemical complexities of the molecular 
gas in galaxies. The 3 mm wavelength band is equally 
well suited, and past spectral surveys have provided im- 
portant information on the chemical and physical prop- 
erties of Giant Molecular Clouds (GMCs) in the Milky 
Way (iJohansson et al.iri984 119851: [Cummins et"aHll986l: 
lTurnerl ll989L One major disadvantage of most spectral 
line surveys is that the data are assembled from many ob- 
servations with varying pointing accuracy and with po- 
tentially systematic calibration problems. The observa- 
tions described here have been obtained simultaneously 
over the full spectral band for each galaxy, and hence 
many systematic problems are eliminated. In this paper 
we present the first 3 mm spectral scans of the central 
regions of 10 galaxies. 
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2. OBSERVATIONS AND DATA REDUCTION 

The observations reported in this paper were ob- 
tained in spring 2008 using the Redshift Search Re- 
ceiver (RSR) on the Five College Radio Astronomy Ob- 
servatory (FCRAO) 14-m telescope. The RSR is a 
sensitive, ultra-broad bandwidth receiver/spectrometer 
(jErickson et al.1 I2007T ) developed at the University of 
Massachusetts as a facility instru ment for the 50 -m di- 
ameter Large Millimeter Telescope (jSchloe rb 2008) . This 
instrument was designed primarily to measure the red- 
shift of distant, dust-obscured galaxies. The RSR is a 
dual polarization and dual beam instrument. The four 
broadband receivers cover instantaneously the frequency 
range 74-111 GHz. A high speed Faraday rotation beam 
switch, operating at 1 kilohertz, is used to overcome 
the 1// noise originating in the front-end monolithic mi- 
crowave integrated circuit (MMIC) amplifiers. Following 
the MMIC amplifiers, two wideband mixers convert each 
receiver band to two intermediate frequency (IF) chan- 
nels. After further conversion and amplification, the IF 
signal is passed into an analog auto-correlation spectrom- 
eter. Each analog correlator has a bandwidth of 6.5 GHz 
and there are six correlators for each receiver polariza- 
tion. To obtain the best frequency resolution, we do not 
apodize the lag domain signal before transforming to the 
frequency domain. Without this apodization, as with 
other auto-correlation based spectrometers, a ringing ef- 
fect can be seen in the baselines around strong, narrow 
lines. The RSR has been designed to detect weak rel- 
atively broad lines, in which case this ringing is not a 
problem. In several of the spectra presented here, the 
13 CO line is sufficiently strong, that this ringing adds 
additional baseline noise near the line. The RSR has an 
instantaneous bandwidth of 37 GHz with a resolution of 
31 MHz, or a velocity resolution of approximately 100 
km s _1 in the 3 mm wavelength band. 

The RSR was commissioned on the FCRAO 14-m tele- 
scope in 2007 and 2008 and used for several initial sci- 
ence projects. During the time of the observations re- 
ported here, 12 of the final 24 spectrometers were avail- 
able, which permitted beam switching with the full 37 
GHz bandwidth in a single polarization. During commis- 
sioning observations, there was a small hardware issue 
(which has since been diagnosed and fixed) which pro- 
duced anomalous noise at approximately 92 GHz, and a 
small region around this frequency has been blanked in 
many of the spectra presented here. 

Observations were obtained at the positions of the 10 
galaxies listed in Table 1. Distances in Table 1 for the 
nearby galaxies, NGC 253, Maffei 2, IC342 and M82, 
are from IKarachentsevI (|2005D and distances for the more 
distant galaxies are from the NASA/IPAC Extragalactic 
Database. The galaxies were selected primarily because 
they were well studied and had relatively bright molec- 
ular emission lines, important for verifying the perfor- 
mance of the RSR. Since our observations cover the en- 
tire 3 mm wavelength window with uniform sensitivity, 
we can study the previously detected lines and all other 
lines within frequency range of the RSR. We added to 
the list of bright line galaxies some additional weak-line 
galaxies to further test the performance of the RSR and 
many of these additional galaxies are known to host an 
active galactic nucleus ( AGN) . Also included in the table 



is the total integration time spent on each galaxy. The 
observations were taken over varying atmospheric con- 
ditions, so our broad aims were to achieve a relatively 
uniform sensitivity of about 1 mK or to achieve a signal 
to noise in the 13 CO line of 10. NGC 3079 and NGC 
6240 did not show strong lines, so we spent additional 
time on these sources to lower the noise to better than 
0.5 mK. The observations for each galaxy were obtained 
over several observing sessions and combined together. 
The pointing and calibration was repeatedly checked by 
observations of the continuum emission from planets and 
quasars. Pointing offsets were never more than a few arc- 
seconds, a small fraction of the beam size, and the overall 
flux calibration is repeatable to be better than 10 %. 

Also included in Table 1 are some brief notes regarding 
the properties of the central regions of the galaxies in our 
sample. The molecular emission in many of these galax- 
ies is dominated by their nuclear starburst. Both NGC 
6240 and Arp 220 are ultra luminous infrared galaxies 
(ULIRGS). NGC 1068, NGC 4258 and NGC 6240 all 
have Seyfcrt 2 nuclei; however, the emission we observe 
in NGC 1068 is likel y dominated by the sur rounding nu- 
clear starburst ring (jSchinnerer et alJl2000T ). NGC 4258 
is a weak AGN, but has a pair of ra dio jets that may b e 
influencing the molecular emission (|Krause et alJ feOO?). 
NGC 6240 is the merger of two Seyfert 2 host galax- 
ies where the two AGNs are separated by less than 2" 
and the molecular emission is concen trated in a sma ll 
region centered on the Seyfert nuclei (jlono etaD[2Q03). 
Finally we note that NGC 3690 is in the process of merg- 
ing with IC 694 (the system is called Arp 299 or Mrk 
171). The nuclei of NGC 3690 and IC 694 are separated 
by only g 20 arcseconds a nd both nuclei have an AGN 
(jGarcia-Marin et al.ll200(l . The 13 CO and HCN emis- 
sion from the IC 694 nucleus is much stron ger th an that 
from the NGC 3690 nucleus (jAalto et alJll997t ): thus, 
although we are centered on NGC 3690, the emission 
may have a contribution from IC 694. Unfortunately, 
the emission from both nuclei is at approximately the 
same velocity, so velocity cannot be used to separate the 
emissions within our telescope beam. 

One limitation of the RSR for the observations pre- 
sented in this paper is the relatively small throw of the 
beam switch. The reference beam is only offset by 4.34 
arcmin in azimuth. Thus, for the largest galaxies, such 
as IC 342 or NGC 253, the reference beam may con- 
tain weak emission from molecular clouds in the galactic 
disk. However, the 12 CO emission falls off sharply from 
the central regions, and even in the largest galaxies, the 
emission in the reference beam is more t han 10 times 
weak er than that at the center of the galaxy ([Young et al.l 
1995). As a check on line fluxes we compared the inte- 
grated intensity of the 13 CO emission, which of the lines 
we observe has the greatest potential for disk contam- 
ination, in the four galaxies (IC 342, NGC 253, NGC 
1068 and NGC 3079) that we re in common w i th the 
study of iPaglione et al.l (|2001h . iPaglione et all (|2001l ) 
also used the FCRAO 14 m telescope; however they ob- 
tained their observations by position switching well off 
the galaxy. We find in all five galaxies that the integrated 
intensity we measure is in good agreement with that of 
IPaglione et al.l (|2001l k thus, we believe that the line inte- 
grated intensities that we report have not been strongly 
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affected by our beam switching observation mode. 

The data have been analyzed using software developed 
by G. Narayanan specifically for RSR data reduction, 
and this software is a complete data reduction package 
with interactive graphics implemented using the Python 
programming language. For the data for each galaxy, 
linear baselines were subtracted from each of six spec- 
tral segments independently. The six segments were then 
combined together to form a spectrum covering 74 - 111 
GHz. The complete spectrum of each of the 10 galaxies 
is presented in Figure 1. Note that due to the problems 
mentioned earlier, in many of the galaxies a small spec- 
tral window around 92 GHz has no usable data. 

The RSR does not have Doppler tracking ability; thus, 
the spectra have a frequency axis in the local rest frame 
of the telescope. Due to changes in the Earth's motion 
relative to the sources, observations made at different 
times and on different days have slightly different fre- 
quency offsets, but these offsets arc less than 5 MHz, 
much smaller than the 31 MHz resolution of the spec- 
trometer. The intensity scales of all spectra shown are 
T* A (antenna temperature corrected only for atmospheric 
losses). The main beam efficiency of the 14-m telescope 
using the RSR has been measured to vary from 0.63 at 89 
GHz to 0.51 at 107 GHz. The decline in efficiency at high 
frequencies is due to the increased amount of power lost 
to the error pattern of the telescope. The beam size also 
varied with frequency, with a HPBW of approximately 
70" at 74 GHz decreasing to 47" at 111 GHz. 

Since the galaxies observed vary enormously in dis- 
tance, the regions probed by our RSR observations vary 
from spatial scales of approximately 1 kpc for the nearest 
galaxy to 18 kpc for the most distant galaxy. However, in 
many of these galaxies the molecu lar emission is centrally 
concentrated (|Young et al.1119951 ). so that the observed 
emission is likely dominated by gas in the central ~1 
Kpc region. It is worthwhile noting that the beam size 
variations over the RSR bandpass could distort line ra- 
tios if the emission were extended. However, if the bulk 
of the emission arises from a nearly unresolved source, 
then the line ratios will not be affected. 

3. RESULTS 

3.1. Line-integrated Intensities 

The spectra shown in Figure 1 display a number of 
strong lines that are in common to almost all of the galax- 
ies observed. In addition to the strong lines of 13 CO, 
HCN, HCO+, HNC and CS, many additional much 
weaker lines are also detected. T he spectral scans of 
Orion KL and Sgr B2 (|Turnerl ll989) were used as a guide 
to compile a list of spectral lines and line blends that 
are likely to be detectable in these galaxies and from this 
compilation we have identified all of the spectral features 
in the 10 galaxies. Lines from CO, CS, HCN, HCO+, 
HNC, N 2 H+, C 2 H, HNCO, CH 3 OH, HC 3 N, C 3 H 3 , SO 
and CH3C2H are detected in one or more of these galax- 
ies. In addition, at least one of the hydrogen recombi- 
nation lines was detected in M82, NGC 1068 and NGC 
4258. 

One of the most important advantages of obtaining 
each 3 mm wavelength spectrum simultaneously with the 
RSR is the perfect pointing registration in all spectral 
lines. These galaxies all have complex small-scale struc- 



ture, so even small pointing errors can seriously affect 
line ratios. In addition, because the lines were obtained 
simultaneously with a single instrument, the relative cal- 
ibration of the spectral line intensities is very good. 

The intensities of the spectral features were determined 
by fitting each spectrum with a template of spectral lines 
assuming a common line width and redshift. The list 
of lines was developed from the stro ng lines seen in the 
spectral scans of Orion and Sgr B2 (j Turner! [19891 . and 
include all plausible lines, even if the chemistry is differ- 
ent than these Galactic sources. The benefit of such a 
simultaneous fit is that the line width and redshift can be 
well constrained by the stronger spectral lines, permit- 
ting more accurate determination of the intensities of the 
weaker lines. Of course this assumes that the emission 
from the various molecular species are distributed simi- 
larly so they share a common line width and line velocity. 
Examination of the residuals to the template fits finds 
them consistent with the baseline noise, thus, we believe 
that these assumptions are satisfied. The 33 spectral fea- 
tures that were fit are listed in Table 2, which include the 
isotopic lines of H 13 CN, H 13 CO+, HN 13 C, and C 34 S. Al- 
though these isotopic lines were not detected in any of 
the galaxies, they were included in the simultaneous fits 
to provide upper limits on the isotopic emission. Since 
hydrogen recombination lines were possibly detected in 
three of these galaxies, the entire series of hydrogen re- 
combination lines were included in the fits. In several 
cases a single spectral feature was fit to a blend of lines 
(marked blend in Table 2) when the line separation in 
frequency was much smaller than the 31 MHz spectral 
resolution of the RSR. In most cases these were blends 
of lines from the same molecular species. 

Fits to the 33 spectra features are shown superimposed 
on the spectra in Figure 1. All features with any signifi- 
cance are well fit by this line template. Features detected 
at a significance of 3a or greater are labeled in Figure 1. 
The results of our template fitting procedure for all of 
the galaxies are summarized in Table 3. For features 
detected at a significance of 3a or greater, the table pro- 
vides the peak intensity and the la uncertainty on these 
intensities (in parentheses). For features not detected at 
the 3a significance level, the table provides the 3a up- 
per limit to the intensity of the spectral transition or 
blend. Only one of the six hydrogen recombination lines 
was detected in M82, however, two other recombination 
lines were just slightly below the detection threshold and 
remaining two recombination lines were present at the 
2a level. Thus, the presence of hydrogen recombination 
lines is probably secure. In NGC 1068 and NGC 4258 
only the H41a recombination line was detected and in 
NGC 4258 this was only one of the three lines detected. 
In both of these galaxies there is no evidence for any of 
the other recombination lines at comparable strengths as 
the H41a line as would be expected, therefore, the iden- 
tification of this feature at 92.035 GHz as a hydrogen 
recombinations line is suspect. The only other plausible 
identification for this feature is a blend of CH 3 CN lines 
from the J = 5-4 transition, however if correct, CH 3 CN 
would be surprisingly strong in these two galaxies. 

The only spectral line detected in all 10 galaxies is 
HCO+. 13 CO was detected in nine galaxies, HCN, HNC 
and CS were detected in seven galaxies and C 18 was 
detected in six galaxies. Because our sample of galaxies 
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vary in distance from 2.9 to 101 Mpc, the absolute de- 
tection threshold is not the same for these galaxies, how- 
ever, in all subsequent analysis we utilize line ratios. The 
richest spectrum is that for NGC 253, where 16 of the 33 
spectral features listed in Table 2 were detected. We note 
that there have been many spectroscopic observations of 
these same galaxies, and a number of the spectral lines 
in Table 2, not detected in our survey, were detected in 
more sensitive measurements focused on single spectral 
lines. In only six galaxies, Arp 220, NGC 1068, NGC 253, 
M82, IC 342, and Maffei 2, five or more spectral features 
are detected. In NGC 6240 the only line detectable is 
HCO+, and in two other galaxies, NGC 3690 and NGC 
4258, we detect HCO+ but not HCN. In the following 
section the template fits are used to examine several line 
intensity ratios and investigate variations in these ratios 
from galaxy to galaxy. 

3.2. Measured Line Ratios 
3.2.1. Isotopic Ratios 

We first examine the CO isotopic ratios. In the six 
galaxies where both 13 CO and C 18 were detected, the 
average ratio 13 CO/C 18 was 4.5. The spread in this 
ratio was relatively small, ranging from 3.3 in NGC 
3079 to 7.8 in Maffei 2. These measured ratios are 
similar to those fo und for massive G MCs in the disk 
of the Milky Way dPenzias et al.lfl972f) and for Sgr B2 
(jCummins etal! 11986ft . A~ surve y of luminou s star- 
forming regions in the Milky Way (jLoenenl 12001 ) found 
13 CO/C 18 ratios between 3 and 16, with the average 
region having a ratio of about 6, again very similar to 
the values measured in our galaxy sample. This ratio 
can be affected by optical depth; however, this lowest 
transition of CO is intrinsically weak, so it is expected 
that these isotopic lines are optically thin. Thus, this 
ratio likely reflects the abundance ratio of the isotopes 
of atomic carbon and oxygen. In the Milky Way both 
the 12 C/ 13 C ratio and the 16 Q/ 18 Q ratio hav e a strong 
dependence on galactic radius (| Wilson! |1 9991). both in- 
creasi ng with increasing galactic radius. iWilson fc Rood! 
(1994) suggest that both 13 C and 18 are secondary nu- 
clear products, while 1 C and 16 are primary nuclear 
products. Ratios of primary to secondary products, such 
as 12 C/ 13 C, are expected to be a good indicator of the 
chemical evolution of the Galaxy. However ratios of sec- 
ondary products, such as 13 CO/C 18 0, should not vary 
greatly as is observed. 

We did not detect the J=l-0 transitions of H 13 CN, 
H 13 CO+, or HN 13 C or the J=2-l transition of C 34 S in 
any of the galaxies in our sample, however, these limits 
provide information on the optical depth of the main 
lines. The most sensitive limits on the isotopic line 
strengths were set in NGC 253 and M82. In NGC 253 
the 3-<T lower limits on the isotopic intensity ratios were 
HCN/H 13 CN > 19, HCO+/H 13 CO+ > 16, and CS/C 34 S 
> 8. In M82 the 3-cr lower limits were HCN/H 13 CN > 12, 
HCO+/H 13 CO+ > 19, and CS/C 34 S > 5. If we assume 
iso topic abundan ces measured or NGC 253 summarized 
bv lOmontl (120071) . we conclude that the main isotopic 
lines of HCN, HCO+, HNC and CS have optical depths 
< 3. 

3.2.2. HCN/ 13 CO Ratio 



Both CO and HCN emission are often used as tracers 
of the molecular gas in galaxies (|Omontll2007l ). Since the 
transitions of CO and HCN have quite disparate crit- 
ical densities, they may be expected to trace different 
density molecular gas, and thus their ratio is a measure 
of the dense gas fraction ([Paglione et al.l [T998h . 12 CO 
is outside the frequency range of the RSR for these very 
low redshift galaxies; however, we believe that 13 CO may 
be a much better probe of the gas. The 13 CO emission 
likely arises from higher column density molecular gas 
than CO, and much of the 13 CO emission may be pro- 
duced in the same gas responsible for the HCN emission. 
Thus, the HCN/ 13 CO ratio may be a good diagnostic of 
the density of the molecular gas most directly connected 
to the star formation activity in the nuclear starburst 
regions of these galaxies. We find that the HCN/ 13 CO 
ratio varies significantly from galaxy to galaxy, with the 
largest ratio of 2.2 found in Arp 220 and the smallest 
ratio found in NGC 4258, where the 3-cr upper limit was 
only 0.30. We discuss later how this ratio can provide an 
important constraint on the mean gas density. 

3.2.3. HCO+/HCN Ratio 

Early models of the chemistry of X-ray irradiated 
gas suggested that HCO + woul d be under-abundant in 
gas near a hard X-r ay source (jLepp fc Dalgarncl 119961 : 
iMalonev et al.l [19961) . These results led to the sugges- 
tion that the HCO + /HCN ratio may distinguish be- 
tween PDR and XDR dominated gas and thus between 
central regions of galaxies dominated by either star- 
bursts or AGNs. Subseque nt obse rvations (IKohno et al.l 
20011: IGrada-Carpio et ail 120061: llmanishi et al.l 120071: 
Baan et al.l |2008|) seemed to confirm this prediction. 
These observations also showed an anti-correlation be- 
tween the observed HCN/CO ratio and the HCO+ /HCN 
ratio. However, th e mor e rece nt theoretical work o f 
IMeiierink fc Spaansl (<2005D and IMeiierink et al.l (120071) 
places this interpretation in question. IMeiierink et al.l 
(2007) showed that in the high density gas, the 
HCO + /HCN ratio is larger in XDR regions than in 
PDR regions, just the opposite of that suggested by 
Lcpp & Dalgarno (1996). They also noted a strong den- 
sity dependence for this ratio in both PDR and XDR 
gas. 

In Figure 2 we plot the HCO+/HCN ratio versus 
the HCN/ 13 CO ratio. We find a similar trend to 
that found in previous observational papers, that is, 
the HCO + /HCN ratio is inversely correlated with the 
HCN/ 13 CO ratio. The agreement of our result with past 
papers is not surprising, since there is considerable over- 
lap between our galaxy sample and those in the previous 
studies. We will address the role density plays in these 
ratio trends later, however, NGC 3690, NGC 4258 and 
NGC 6240 are quite anomalous in having unusually large 
ratios of HCO + /HCN, and variations in chemistry are 
almost certainly needed to explain these ratios. NGC 
6240 is not plotted in Figure 2 because we did not de- 
tect either HCN or 13 CO, but the ratio of HCO+/HCN 
is greater than 2.3 (3 sigm a limit). NGC 3690 and NGC 
6240 were both part of the Uuneau et al.1 ()2009l ) study of 
ULIRGS, and these galaxies were among those with the 
largest HCO+ /HCN ratio. 
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3.2.4. HNC/HCN Ratio 

The HNC/HCN ratio is a tracer of gas temperature 
(jAaltol [2008T ) . It is well known that in cold clouds the 
colum n density ratio of HNC /HCN can be as large as 3 
to 10 ([Churchwell et al.|[l98l . while in warm GMCs, the 
column density ratio is less tha n 1, and can be as small 
as 0.015 (Goldsmit h et al. l ll981f ). The strong dependence 
of this ratio on temp erature is a predic tion of molec- 
ular cloud chemistry (Schilkc et al. 1992). Since these 
lines are similarly excited and have similar yl-coefficients, 
when they are optically thin the observed emission ra- 
tio is proportional to the column density ratio. For the 
galaxies in our sample with detections in both lines, the 
observed ratio of HNC/HCN is confined to a narrow 
range of values between 0.43 and 0.69, n a rrowe r than 
that the galaxies observed by lAalto et al.1 ([2002). The 
ratios are consistent with the emission arising from warm 
molecular gas. 

3.3. Comparison of M82 and NGC 253 

The two galaxies with the strongest emission lines are 
M82 and NGC 253 which enable a comparison of emis- 
sion from other molecular species, such as C2H, N9H+. 
HNCO, CH3OH, HC 3 N, C3H2, SO and CH 3 C 2 H, that 
had only a limited number of detections in our galaxy 
sample as a whole. In both of these galaxies the emis- 
sions observed are likely dominated by molecular clouds 
associated with their nuclear starbu rsts. However there 
are differences, lAladro et al.l ()2011[ ) suggested that the 
starburst in M82 is more evolved and the chemistry is 
dominated by PDR regions, while the chemistry in the 
younger starburst, NGC 253, is dominated by shocks. 

Notable differences between these two galaxies are 
the detection of HNCO, CH 3 OH, N 2 H+ and HC 3 N in 
NGC 253 but not in M82, although these lines would 
have been readily detected if they had the same line 
strength relative to HCN as they have in NGC 253. Both 
HNCO and CH3OH are thought to have a shock ori- 
gin, as t hese molecules are be lieved to form on grain 
surfaces (jTideswell et al.l 120101 ) an d then ejected from 
the grains by low velocity shocks ([Martin et al.l [20081. 
Thus, the presence of these lines in NGC 253 agrees 
with the idea that shocks are an important component 
to the c hemis try in the nuclear regions of this galaxy. 
lOmontl (|2007f ) also noted that SiO and CH 3 CN, both 
thought to have a shock origin, are more abundant in 
NGC 253 than M82. We note that HNCO and CH 3 OH 
are even stronger relative to HCN in IC 342 and Maffei 2, 
both whose chemist ry may also be dominated by shocks 
(|Aladro et al.1120111 ) . A similar conclusion c oncerning IC 
342 was reached bv lMeier fe Turned (|2005h . The largest 
CH3OH/HCN ratio was measured in NGC 3079, and 
thus the chemistry in NGC 253, IC 342, Maffei 2 and 
NGC 3079 may all be strongly influenced by shocks. 
HC 3 N was also detected in NGC 253 and not in M82. 
The only other galaxy where HC3N was detected was 
Arp 220, where the ratio of HC3N/HCN was about three 
times larger than in NGC 253. Base d on the stron g emis - 
sion in Arp 220 ([Aalto et alJl2002f ). lAalto et al.l (|2007f l 
speculated that strong HC3N emission might be indica- 
tive of young starbursts and this perhaps may explain 
the stronger emission in NGC 253 relative to the older 
starburst M82. 



The other notable difference is the detection of 
CH3C2H in M82 but not in NGC 253, and the 
CH3C2H/HCN ratio measured in M82 is three times 
larger than the 3cr upper limit established in NGC 253. 
In fact, other hydrocarbon molecules, such as C2H, 
C3H2 and CH3C2H, have strengths relative to HCN 
much larger in M82 than in NGC 253. The enhance- 
m ent of these hydroca rbon molecules in M82 is discussed 
in iFuente et all ([2005D and they suggested that photon- 
dominated chemistry was important for the production 
of hydrocarbons in the central region of M82. Thus the 
pattern of lin es seen in M82 and NGC 253 supports the 
suggestion of lAladro et al.l ([201 ll ) that the chemistry in 
NGC 253 is influenced by shocks, while the chemistry in 
M82 is more PDR dominated. 

3.4. Starburst Versus AGNs 

From our sample of galaxies, we can also examine 
the differences between the starburst and AGN galax- 
ies. We have already noted that NGC 3690, NGC 4258 
and NGC 6240, all with AGNs, have anomalously large 
HCO+/HCN ratio. We have formed composite spectra 
of a sample of starburst galaxies (NGC 253, Maffei 2, IC 
342, M82 and Arp 220) and galaxies with AGNs (NGC 
1068, NGC 3690, NGC 4258 and NGC 6240). The emis- 
sion in each spectra was first normalized to the strength 
of the HCO + line (the only line seen in all the galaxies), 
then the frequency axis was adjusted by a factor of 1+z, 
and finally the spectra were averaged using equal weights. 
Equal weighting was used so that no galaxy dominated 
the composite spectrum, although this meant that the 
weaker line galaxies added more noise than the strong 
line galaxies to the composite spectra. The composite 
spectra were fit using the same procedure described ear- 
lier for the individual galaxies and, as before, a 3a cutoff 
was used to define line detections. Unfortunately, the 
AGN composite consisted primarily of weak-line galax- 
ies, so fewer lines were detected in the composite spec- 
trum relative to the starburst composite. The most in- 
teresting ratio is the HCO + /HCN ratio and in the star- 
burst galaxies this ratio has a value of 0.71 ±0.06, while 
in the AGN galaxies, the ratio was 1.9±0.3, significantly 
larger. We can also compare the HCN/ 13 CO ratio, which 
for the starbursts was 1.11±0.07 and for the AGNs was 
0.64±0.09. Ratios of HNC/HCN and CS/HCN were sim- 
ilar in the two samples. 

In the AGN galaxies, NGC 1068 has the strongest lines. 
Earlier we noted that none of the observed line ratios 
distinguishes NGC 1068 from the sample of starburst 
galaxies. The emission observed in NGC 1068 may be 
dominated by the starburst ring and not gas closely as- 
sociated with the Seyfert nucleus. Removing NGC 1068, 
we can form a new AGN composite spectrum from the 
remaining three AGN galaxies. Although HCN was not 
detected in any of the individual galaxies, it is readily 
detected in their composite spectrum. In this new com- 
posite AGN spectrum, we find ratios of HCO+/HCN = 
3.4±0.9 and HCN/ 13 CO = 0.45±0.12. Not surprising, 
these ratios are even more extreme relative to the star- 
burst sample than the AGN composite including NGC 
1068. 

4. MOLECULAR GAS MODELS 



6 



Snell et al. 



4.1. Modeled Line Ratios 

The line ratios discussed in the previous section de- 
pend on both the physical (density, column density, and 
temperature) and chemical (molecular abundance) prop- 
erties of the gas. The information available in our spectra 
is inadequate to constrain all the needed gas properties. 
We have consequently modeled the line ratios with a very 
simplistic approach in which we assume that the emis- 
sion arises from a collection of identical molecular cloud 
cores which are characterized by a single temperature 
and density, with equal abundances of the various molec- 
ular species. Molecular abundances were taken from 
the well studie d Galactic GMC co res M17 and Ceph A 
summarized bv iBergin et al.l (|1997f ). We further assume 
abundance ratios of CO/H 2 = lxlO" 4 and 13 CO/CO = 
2xl0 -2 ; the resulting molecular abundances relative to 
H 2 are given in Table 4. An E/A ratio of 1 for CH 3 OH 
and an ortho/para ratio of 3 for C3H2 were also assumed. 
The assumed m olecular abundances are very similar to 
those derived bv lMartm et al.1 (|2006f ) for NGC 253 based 
on their 2 mm wavelength spectral survey (see compari- 
son in Table 4). 

The emission from an ensemble of GMC cores was 
mode led with the RADEX code dVan der Tak et al.l 
[2001 using the Einstein ^-coefficients and collision 
rate coeffi cients from the Leid en Atomic and Molecular 
Database (Schoicr et al. 2005). To compute the line in- 
tensities, besides the molecular abundances, the gas den- 
sity, n, gas temperature, T, and the H2 column density 
per unit line width, N(H2)/Ai>, need to be specified. Our 
standard GMC core model assumes T = 35 K, n = 1 x 10 5 
cm" 3 , and N(H 2 )/Av = lxlO 22 cm" 2 (km s -1 ) -1 . Line 
width s for GMCs in t he Galactic disk are around 3 km 
s _1 (jSnell et al.l 119841 ) . leading to a H 2 column density 
of 3xl0 22 cm -2 . Molecular clouds near the center of our 
Gala xy, such as Sgr B2, ha ve line widths around 20 km 
s _1 (jCummins et al.l 119860 . and for this line width the 
modeled H2 column density is 2xl0 23 cm~ 2 . 

The line ratios computed from our standard model 
are compared with the observations of NGC 253, the 
galaxy with the richest spectrum, in Figure 3. Since 
HCN usually has the strongest emission of any of the 
'high dipole moment' molecules, all intensity ratios are 
computed relative to HCN. Multiple lines were observed 
for several molecular species and for the line ratio analy- 
sis the 87.9252 GHz line for HNCO, the 96.7414 GHz line 
for CH3OH and the 99.2999 GHz line for SO were used 
in forming the intensity ratios. For HC3N, the two lines 
at 81.8815 and 100.0764 GHz, which were in the least 
line confused regions of the spectrum, were summed to- 
gether. The line ratios for our standard core model agree 
well with those observed for NGC 253. The largest dis- 
crepancies between the observed and modeled line ratios 
is less than a factor of two. The observed ratios in this 
galaxy are well fit using this simple model. 

Our model can also be used to address the sensitivity of 
the various line ratios to density and temperature. The 
line ratios were computed from the standard model for 
densities of 1 x 10 4 and 1 x 10 6 cm -3 and the resulting line 
ratios shown in Figure 3. For our assumed abundances, 
a density of 1x10 s cm~ 3 is the better fit to the observed 
line ratios in NGC 253. Figure 3 illustrates that several 
of the line ratios, due to differences in the critical density 



of the lines forming the ratio, are highly sensitive to den- 
sity. As expected, one of the values most affected is the 
HCN/ 13 CO ratio; however, in addition, the line ratios 
CH3OH/HCN and HNCO/HCN are also density sensi- 
tive. The temperature sensitivity of these ratios can also 
be explored, and in Figure 4 the observed line ratios in 
NGC 253 are compared to those derived from our stan- 
dard model with density 1x10 s cm~ 3 and with varying 
temperatures. Temperature has a much smaller effect 
on the line ratios than does density, largely because the 
molecular transitions we observed all arise from similar 
energy levels above the ground rotational state. 

NGC 253 has been widely studied and there are a 
number of estimates of the density of the molecular gas. 
Ba sed on fitting multip le li nes of CS, HCN and HCO+ 
led lMartin et al.1 (|2005h and lKnudsen et al.l (|2007t ) to es- 
timate a density of approximately 2xl0 5 cm -3 , higher 
than an earlier estimat e based o n H2CO emission 
(jHuttemeister et al.lll997h . Recently. iBavet et al.1 (|2009[ ) 
fit the CS emission from NGC 253 with a two-component 
model consisting of a cool, lower density component and 
a warm, higher density component. For the two veloc- 
ity features in NGC 253, they found densities for the 
cooler component of order 2xl0 4 and for the warmer 
component of order 2xl0 6 cm~ 3 . The density in our 
st andard model agrees well with the single densit y fits 
of iMartm et all (|2005l ) and iKnudsen et all (|2007t ) and 
is consistent with the average of the components fit by 
iBavet et alj ([2009). Based on this density, the abun- 
dances we assumed in our standard model, including 
those molecules (HNCO and CH 3 OH) that may be pro- 
duced in shocks, are consistent with the observations in 
NGC 253. A better determination of the physical prop- 
erties of the molecular gas is needed to refine the abun- 
dance estimates. 

Our model can be used to examine the anti-correlation 
between the HCO+/HCN and 13 CO/HCN ratios shown 
in Figure 2. Overlaid on the data in Figure 2 is shown 
a line connecting the results for five models in which 
only the density is varied, from lxlO 4 to lxlO 6 cm~ 3 . 
The comparison between data and models suggests that 
much of the variation observed in the 13 CO/HCN ratio 
can be explained if the average density of the molecu- 
lar gas varies among these galaxies. For densities below 
lxlO 4 cm -3 , it is difficult to produce detectable HCN or 
13 CO emission without having unrealistically large col- 
umn densities of these molecular species. Although the 
HCO + /HCN ratio is not strongly affected by density, 
some of the anti-correlation found between this ratio and 
the 13 CO/HCN rati o can be expla i ned b y variations in 
the average density. Uuneau eta l. (2009) came to simi- 
lar conclusion and the sensitivity of these line ratios to 
density is also appare nt in the modeling presented by 
iMeiierink et al.1 (|2007f ). Based on the line ratios, and ig- 
noring any chemical differences, the data suggest that the 
molecular gas in Arp 220 is the densest and the molec- 
ular gas in NGC 3079 the least dense. However, density 
alone cannot explain the large HCO + /HCN line ratios 
measured in NGC 3690, NGC 4258 and NGC 6240, a ll 
with AGNs. In the models of IMeiierink et al.l (|2007fl . 
only the XDR models predict HCO + /HCN ratios greater 
than one, thus, if these XDR models are correct, it is 
possible that the AGN has an influence on the chemistry 
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and is responsible for the enhanced HCO + abundance 
in the central regions of these three galax ies. However, 
the la rge sample of ULIRGS presented by Uuneau et al.l 
(2009) show no correlation between the HCO + /HCN ra- 
tio and whether these galaxies are believed to be AGN 
or star formation dominated, putting some doubt in this 
AGN hypothesis for the enhanced HCO + emission. 

The comparison between models and observations can 
be extended to the line ratios measured for IC342, M82, 
Arp 220, Maffei 2, NGC 1068 and NGC 3079. These com- 
parisons are shown in Figure 4. Since density has the 
greatest influence on the line ratios, only models with 
varying density at a fixed temperature of 35 K are in- 
cluded in this figure. Based on model comparisons shown 
in Figures 2 and 4, IC 342, Maffei 2 and NGC 3079 would 
be expected to have the lowest average gas density of 
~10 4 cm -3 . This low average density is borne out by 
the other density sensitive ratios. With this density, we 
find the abundances from our standard core model repro- 
duce well the measured line ratios in these three galaxies. 
Similarly, the average density for M82 is expected to be 
somewhat larger (^3x10 cm -3 ) and assuming this den- 
sity we confirm that the abundances of HC3N, CH3OH, 
N 2 H + , and HNCO relative to HCN are smaller than in 
our standard model and for those found in NGC 253. 
Again based on the result shown in Figures 2 and 4, NGC 
1068 is expected to have a somewhat larger average den- 
sity (~10 5 cm -3 ), while Arp 220 is expected to have the 
highest average density (~10 6 cm" 3 ). As suggested ear- 
lier, we find an overabundance of HC3N relative to HCN 
in Arp 220 for any of the modeled densities. Thus if we 
allow for the average molecular gas density to vary from 
galaxy to galaxy, our standard model with fixed abun- 
dances provides a reasonable good fit to the measured 
ratios in these galaxies. 

There have been a number of estimates of the av- 
erage molecular gas density for galaxies in our sam- 
ple. Most of these estimates utilize multiple transi- 
tions of a single molecular species, thus they are in- 
dependent of derived abundances. However, despite 
this, estimates for a single galaxy can still differ by 

order of magnitude . De nsity estimates for I C 34 2 

(|Hiittemeister et al.l 119971: IMeiierink fc Spaansl I2005D. 
Maffe i 2 (IHiittemeister et al.l 119971). M8 2 ()Navlor et al l 
I2010D . and Arp 220 (jGreve et"alTl2009D are in reason- 
able agre ement with our e stima tes. The HCN line SEDS 
shown in iKnudsen et all (|2007f) are also consistent with 
the trend in density variations for the galaxies M82, NGC 
253 and Arp 220. Fitting the gas in the central regions 
of these gal axies with a single dens ity is an over s impli - 
fication and lBavet et~aTl (|2009[ ) and lAladro et al.l (|2011f l 
have fit multiple-density components for NGC 253, M82, 
IC 342, Maffei 2 and NGC 1068. With the exception 
of M82, their column density weighted average density 
is consistent with our findings. For M82, both find a 
column density weighted d ensity about a n ord er of mag- 
nitude higher than that of Na vlor et al.l (|2010f ) . 

4.2. Molecular Mass and Beam Filling Factor 

The simple model we presented in the previous sec- 
tion, in which the molecular emission observed arises 
from an ensemble of identical cloud cores, can be used to 
estimate the beam filling factor and mass of molecular 
gas, independent of the Xco-factor. Since the distances 



of the galaxies in our sample vary enormously, the re- 
gion probed ranges from just the central starburst core 
to nearly the entire galaxy. The linear diameter of the 
FCRAO beam at the frequency of the 13 CO line for the 
distance of these galaxies (Db) is given in Table 5. 

The beam surface area filling factor (f are a) is related to 
the observed antenna temperature of the line corrected 
for the main beam efficiency ( Tm b ) , the line intensity 
from our model (T mo d e i), the assumed core line width in 
the model (Av mo d e i) , and the observed line width of the 
galaxy (Av oos ), which is dominated by galactic rotation, 
by the following equation: 



farea — 
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The observed 13 CO emission is used to estimate the beam 
filling factor from the above relation. The modeled in- 
tensity ratios, discussed previously, did not depend on 
the line width, however, the total integrated intensity, 
and hence the beam surface filling factor, does. We have 
no information on the line widths of cloud cores in these 
galaxies, so we assume a line width of 10 km s _1 for our 
model cores, which corresponds to an H2 column density 
of 1 x 10 23 cm~ 2 . Although a density of 1 x 10 5 cm~ 3 was 
assumed, the emission from 13 CO has very little den- 
sity dependence as long as n > lxlO 4 cm~ 3 . Based on 
the observed line intensities and line widths presented in 
Table 3, the computed beam area filling factor is summa- 
rized in Table 5. Not surprising, the beam filling factor 
is strongly dependent on galaxy distance. For the nearer 
galaxies, the beam area filling factor can be as large as 
0.26, implying that one-fourth of the FCRAO ~50" beam 
is covered with gas with a total gas column density of 
lxlO 23 cm -2 , corresponding to an A v of approximately 
100 mag. Similar beam area filling factors would have 
been found if we had used other molecular species, such 
as HCN or HCO + . However, the results would have been 
more sensitive to the assumed density of the cloud cores. 

The molecular gas mass can also be estimated from our 
model. The mass is given by: 



M = m H2 N model f area d fli 



(2) 



where mjy 2 is the mass of a molecular hydrogen molecule, 
Nmodel is the gas column density in the core model, d is 
the distance to the galaxy, and Qb is the main beam 
solid angle. As long as the 13 CO emission is optically 
thin, the mass determination is independent of many of 
the assumed core model parameters. However, two im- 
portant parameters in deriving masses are (1) the as- 
sumed core gas temperature, as this affects the partition 
function, and (2) the assumed 13 CO/H 2 ratio of 2xl0~ 6 . 
The mass of molecular hydrogen derived from the above 
equation is summarized in Table 5. Because of the vary- 
ing distance of these galaxies, the mass of the nearby 
galaxies (NGC 253, Maffei 2, IC 342 and M82) includes 
only the gas in the central 1 kpc regions, whereas for 
the most distant galaxies (NGC 3690 and Arp 220) it 
includes the entire molecular mass of the galaxy. 

These mass estimates can be compared with those de- 
rived from CO using an assumed Xco-factor. To make 
certain that we are comparing the mass from comparable 
regions of these galaxies, we have used CO observations 
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also obtained with the FCRAO 14-m teles cope as part 
of th e FCRAO Extragalactic CO Survey (jYoung et al.l 
1995). Observations of all of the galaxies except Maffei 
2 are covered in this survey, and these observations were 
obtained with nearly an identical beam size as our ob- 
servations of 13 CO. The CO fluxes from the survey are 
summarized in Table 5, where the integrated in t ensity 
to flux conversion factor provided in lYoung et al.l (1995) 
was used. 

To estimate the mass from the CO fluxes, an Xqo- 
factor of 1.8xl0 20 cm -2 (K km s -1 ) -1 , was use d which 
is rep resentative of Galactic molecular clouds (jOmont] 
I2007D . The molecular hydrogen gas mass can then be 
written ()Kennev fc Young] [T989) as: 

M{H 2 ) = 7.1 x 10 3 d 2 S CO , (3) 

where d is the distance in Mpc and Sco is CO flux in 
units of Jy km s . The masses derived from this expres- 
sion are presented in Table 5. Mass estimates from CO 
are 1.2 - 4.3 times larger than that derived from 13 CO. It 
is possible to reconcile these masses if we either lower the 
13 CO abundance, increase the gas temperature assumed 
in our core model, or decrease the Xco-factor used. We 
cannot rule out any of these possibilities; however, a more 
likely scenario is that some of the CO emission is arising 
from lower density and lower column density gas that 
does not produce significant 13 CO emission. No matter 
what the reason for the mass difference, we can conclude 
that the mass of dense gas must represent a significant 
fraction of the total molecular gas in these galaxy centers. 

We find that between about 25 and 85 % of the total 
molecular mass in these galaxies must have high density 
and column density. These properties are likely asso- 
ciated with the molecular clouds in the cent ral ~1 kpc 
regions of these galaxies. iBallv et aD (|1987f ) estimated 
the mean density of the molecular gas in the inner 500 
pc radius region of the Milky Way to be 2xl0 4 cm -3 , 
similar to the average density estimated for several of 
the galaxies in our survey. However, the properties of 
the molecular cl ouds in the disk of the Milky Way are 
quite different. iLee et all (|1990t ) mapped a 3 deg2 re- 
gion around I = 24° in both 13 CO and CS. Although, 
no estimate was made of the mean density of the gas, 
the quoted emission ratio of CS/ 13 CO of 0.008 is sub- 
stantially lower than the range of ratios, 0.12 - 1.0, we 
measured for our galaxy sample. The largest CS/ 13 CO 
ratio was observed in Arp 220, which is also modeled to 
have the highest mean molecular gas density. Since this 
ratio is sensitive to the mean gas density, this compari- 
son suggests that a large fraction of the molecular gas in 
the central regions of the galaxies in our sample, and the 
Milky Way, must be dense, while molecular clouds in the 
Milky Way disk have a much smaller fraction of dense 
gas. Studies of i ndividual molecular cl ouds in the disk of 
the Milky Way ()Carpenter et al.lll995l ) confirm this and 
show that high density and high column density cores 
represent only a small fraction of the area and mass of 
the molecular cloud. 

5. SUMMARY 

We have presented the 3 mm wavelength spectra of 
the central regions of ten relatively nearby galaxies. The 
spectra of each galaxy were obtained simultaneously with 



the RSR. The important advantages of obtaining the 3 
mm wavelength spectra simultaneously are the perfect 
pointing registration in all spectral lines and the very 
good relative calibration of the spectral line intensities. 
The only uncertainty in relative line intensities comes 
from the small variations of beam efficiency and beam 
size with frequency. Most of the lines that we detected 
arise from high-dipole moment molecules, which allows 
us to examine the physical and chemical properties of 
the dense gas in these galaxies. Although the galaxies 
in our survey are far from being a homogeneous sample, 
their differences can be exploited to study the extent to 
which various line ratios are found to vary and to what 
extent they require differences in the physical or chemical 
properties of the molecular gas. 

The spectra were fit with a template with 33 spec- 
tral features to determine the line intensities. We de- 
tected 20 of the 33 spectral features in our template, 
and these features arose from 14 different atomic and 
molecular species. Based on the fitted line intensities, 
a number of line intensity ratios were examined. The 
13 CO/C 18 and HNC/HCN line ratios are insensitive to 
density The 13 CO/C 18 ratio was found to have a very 
narrow range of values, from 3.3 to 7.8, supporting the 
idea that both 13 C and 18 are secondary nucleosynthe- 
sis products. The HNC/HCN ratio also was relatively 
constant from galaxy to galaxy, varying between 0.4 - 
0.7 and consistent with the emission arising in relatively 
warm gas. 

The emission from these galaxies was modeled with an 
ensemble of molecular cloud cores, each with a temper- 
ature of 35 K and a molecular hydrogen column den- 
sity per unit line width of lxlO 22 cm" 2 (km s -1 ) -1 . 
The molecular abundances were taken from well-studied 
Galactic clouds, although these ab undances are consis - 
tent with those found in NGC 253 (|Martui et al.l 120061) . 
The sensitivity of various line ratios to density and tem- 
perature was examined and it was found that density 
has the greatest affect. The HCN/ 13 CO ratio is partic- 
ularly density sensitive, and we believe that the large 
variation observed for this ratio may be due to differ- 
ences in the mean gas density in these galaxies. Varying 
the density from 10 4 to 10 6 cm -3 can explain the range 
of HCN/ 13 CO ratios observed and provide consistent re- 
sults with the other density sensitive line intensity ratios. 
Varying the density can also help explain much of the ob- 
served variations in the HCO + /HCN ratio. Thus, a sim- 
ple model with fixed molecular abundances and varying 
density, fits well the 3 mm spectra of most of these galax- 
ies. One of the most notable exceptions is the anoma- 
lously large HCO+/HCN ratios detected in NGC 3690, 
NGC 4258 and NGC 6240, which cannot be fit with- 
out increasing the abundance of HCO + relative to HCN 
relativ e. Based on the XDR model of Mciicrin k et al.l 
420071 ). a plausible explanation is that the AGN activ- 
ity in these galaxies is responsible altering the chemistry 
in these galaxies and producing the large HCO+/HCN 
ratio. 

The cloud core model was used to estimate the filling 
factor of dense, high column density gas in our beam 
and to estimate the mass of dense molecular gas. The 
derived beam filling factors are strongly correlated with 
distance, as our beam varies in linear size from about 
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0.8 Kpc in the nearest galaxies to over 17 Kpc in the 
most distant galaxies. In the nearby galaxies, beam fill- 
ing factors as great as 0.26 are found, suggesting that 
the central regions of these galaxies have a high filling 
factor of dense, high column density gas. We compared 
our mass of dense gas with estimates based on CO and 
an assumed Xco~^ctor and find that the dense gas rep- 
resents approximately 25-85 % of the total molecular gas 
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Fig. 1a. — Three-millimeter spectra for NGC 253 are presented by the binned points. The frequency axis is measured relative to the 
local rest frame of the telescope with no corrections for the Earth's motion. The continuous line (continuous green line in the electronic 
version of the paper) shows the fit of our 33-spcctral feature template discussed in the text. 
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Fig. 1b. — Same as (a) for Maffei 2. 




Fig. lc— Same as (a) for NGC 1068. 
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Fig. Id.— Same as (a) for IC 342. 
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Fig. 1e.— Same as (a) for M82. 
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Fig. If.— Same as (a) for NGC 3079. 
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Fig. 1g — Same as (a) for NGC 3690. 
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Fig. lH.— Same as (a) for NGC 4258. 
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Fig. Ii. — Same as (a) for Arp 220. 
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Fig. Same as (a) for NGC 6240. 
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Fig. 2. — A plot of the observed intensity ratio of HCN/ 13 CO vs. the observed intensity ratio of HCO+/HCN for the galaxies (filled 
circles) where 13 CO was detected. The error bars show the lex uncertainty in the values of these ratios. In the two galaxies where HCN 
was not detected, the open circles with arrows mark the 3cr upper limit to the HCN/ 13 CO ratio and the 3(7 lower limit to the HCO+/HCN 
ratio. The results of our modeling of line ratios (see the text in Section 4.1) are shown by the filled triangles with a line connecting five 
model results with varying density. The models, from left to right, have density of log n = 4.0, 4.5, 5.0, 5.5 and 6.0. 
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Fig. 3a. — The observed intensity ratios (shown as filled circles) for a number of key molecular species is shown for NGC 253. The error 
bars show the lcr uncertainty in the values of the line ratios. Where the molecular species was not detected, the 3<r upper limit on the ratio 
is shown. Accompanying the observed intensity ratios, are ratios computed from our standard core model, see the text; the modeled line 
intensity ratios are connected by a dashed line for log n = 4.0, a solid line for log n = 5.0 and a dotted line for log n = 6.0. 
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Fig. 3b. — Same as (a), except showing models of the line intensity ratios for fixed density of log n = 5.0, and temperatures of 15 K 
(dashed line), 35 K (solid line), and 75 K (dotted line). 
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Fig. 4a. — Same as Figure 3(a) for Maffei 2. 
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Fig. 4b.— Same as Figure 3(a) for NGC 1068. 
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Fig. 4c— Same as Figure 3(a) for IC 342. 
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Fig. 4d.— Same as Figure 3(a) for M82. 
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Fig. 4e.— Same as Figure 3(a) for NGC 3079. 
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TABLE 1 



Galaxy Sample and Observing Times 


Galaxy 


a(2000) 1 


5(2000)! 


Distance 


tint 


Notes 








Mpc 


hr 




NGC 253 


00 47 33 


-25 17 18 


3.9 


3.2 


Starburst 


Maffci 2 


02 41 55 


+59 36 01 


2.8 


1.8 


Starburst 


NGC 1068 


02 42 40 


-00 00 48 


14.4 


3.5 


Starburst, Seyfert 2 


IC 342 


03 46 48 


+68 05 45 


3.3 


8.3 


Starburst 


M82 


09 55 52 


+69 40 47 


3.5 


4.3 


Starburst 


NGC 3079 


10 01 57 


+55 40 46 


19.7 


7.5 


LINER 


NGC 3690 


11 28 32 


+58 33 43 


45 


9.0 


LIRG, AGN 


NGC 4258 


12 18 57 


+47 18 13 


7.9 


3.5 


Seyfert 2 


Arp 220 


15 34 57 


+23 30 09 


76 


6.0 


ULIRG 


NGC 6240 


16 52 58 


+02 24 03 


101 


9.7 


LIRG, Seyfert 2 



1 Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutcs, and arcseconds. 
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TABLE 2 
Fitted Spectral Lines 



Line Number Frequency (GHz) Identification 



1 


74.6446 


H44 a 


2 


79.9127 


H43 a 


3 


81.8815 


HC 3 N 9 - 8 


4 


84.5211 


CH3OH 5(-l) - 4(0) E 


5 


85.3389 


C 3 H 2 2(1,2) - 1(0,1) 


6 


85.4557 


CH3C2H blend 


7 


85.6884 


H42 a 


8 


86.3402 


H 13 CN 1 - 


9 


86.7543 


H 13 CO+ 1 - 


10 


87.0909 


HN 13 C 1 - 


11 


87.3169 


C 2 H blend 


12 


87.4020 


C 2 H blend 


13 


87.9252 


HNCO 4(1,3) - 3(1,2) 


14 


88.6318 


HCN 1-0 


15 


89.1885 


HCO+ 1 - 


16 


90.6635 


HNC 1-0 


17 


90.9790 


HC 3 N 10 - 9 


18 


92.0345 


H41 a 


19 


93.1738 


N 2 H+ 1 - 


20 


96.4129 


C 34 S 1 - 


21 


96.7414 


CH3OH blend 


22 


97.9810 


CS 1 - 


23 


99.0230 


H40 a 


24 


99.2999 


SO 3(2) - 2(1) 


25 


100.0764 


HC3N 11 - 10 (SO 4(5) - 4(4) blended) 


26 


102.5470 


CH 3 C 2 H blend 


27 


106.7374 


H39 a 


28 


108.8939 


CH3OH 0(0) - 1(-1) E 


29 


109.1736 


HC3N 12 - 11 


30 


109.2522 


SO 2(3) - 1(2) 


31 


109.7822 


c 18 o 1 - 


32 


109.9057 


HNCO 5(0,5) - 4(0,4) 


33 


110.2014 


13 CO 1 - 



3 mm Spectra of Galaxies 31 



TABLE 3 

Fitted Line Intensities and Line Widths 



Line Intensities 



Line Number NGC253 


MafTei 2 


NGC1068 


IC342 


M82 


NGC3079 


NGC3690 


NGC4258 


Arp220 


NGC6240 


1 


<2.2 


<2.6 


<1.5 


<1.9 


<1.8 


<0.6 


<0.9 


<1.1 


<0.9 


<0.6 


2 


<2.1 


<2.5 


<1.4 


<1.7 


<1.7 


<0.6 


<0.9 


<1.0 


<0.9 


<0.6 


3 


4 7(0 7) 


<2.5 


<1.4 


<1.7 


<1.7 


<0.6 


<0.9 


<1.0 


<0.9 


<0.6 


4 


<2.0 


<2.4 


<1.4 


<1.7 


<1.6 


<0.6 


<0.8 


<1.0 


<0.9 


<0.6 


5 


2.9(0 7) 


<2.4 


<1.4 


<1.6 


2.7(0.5) 


<0.7 


<0.8 


<1.0 


<0.9 


<0.6 




<2.0 


<2.4 


<1.4 


<1.8 


<1.6 


<0.7 


<0.8 


<1.0 


<0.9 


<0.6 


7 


<2.0 


<2.4 


<1.4 


<1.6 


<1.6 


<0.6 


<0.8 


<1.0 


<0.9 


<0.5 


8 


<2.0 


<2.4 


<1.4 


<1.6 


<1.6 


<0.6 


<0.8 


<1.0 


<0.9 


<0.5 


9 


<2.0 


<2.4 


<1.4 


<1.6 


<1.6 


<0.6 


<0.8 


<1.0 


<0.9 


<0.5 


10 


<2.0 


<2.4 


<1.4 


<1.6 


<1.6 


<0.6 


<0.8 


<1.0 


<0.9 


<0.5 


11 


9.0(0.7) 


<2.4 


2.4(0.5) 


2.0(0.5) 


9.5(0.6) 


<0.8 


<0.9 


<1.0 


<0.9 


<0.7 


12 


4.3(0.7) 


<2.4 


<1.4 


<1.6 


4.6(0.6) 


<0.8 


<0.9 


<1.0 


<0.9 


<0.7 


1.3 


4.6(0.7) 


2.8(0.8) 


<1.4 


2.2(0.6) 


<1.6 


<0.6 


<0.9 


<1.0 


<0.8 


<0.5 


1 1 


37.5(0.7) 


8.2(0.8) 


10.8(0.5) 


11.4(0.6) 


20.1(0.6) 


1.6(0.2) 


<0.8 


<1.0 


2.9(0.3) 


<0.5 


15 


32 6f0 7) 


7.4(0.8) 


8.0(0.5) 


8.7(0.5) 


31.0(0.6) 


1.8(0.2) 


2.5(0.3) 


1.8(0.3) 


1.1(0.3) 


1.2(0.2) 


16 


19.4(0.7) 


5.2(0.8) 


4.7(0.5) 


5.2(0.5) 


8.6(0.5) 


<0.6 


1.2(0.3) 


<1.0 


2.0(0.3) 


<0.5 


17 


3.3(0.7) 


<2.3 


<1.5 


<1.7 


<1.6 


<0.6 


<0.8 


<1.0 


<0.8 


<0.5 


18 


<2.0 


<2.3 


1.5(0.5) 


<1.5 


<1.6 


<0.6 


<0.8 


2.4(0.3) 


<0.8 


<0.5 


19 


5.5(0.6) 


<2.3 


<1.3 


1.7(0.5) 


<1.6 


<0.6 


<0.8 


<1.0 


<0.8 


<0.5 


20 


<1.9 


<2.3 


<1.3 


<1.6 


<1.5 


<0.6 


<0.8 


<0.9 


<0.9 


<0.5 


21 


7.4(0.6) 


3.U(U.8 ) 


<1.3 


4.U(U.0 ) 


<1.5 


U. 6(0.2} 


<0.8 


<0.9 


<0.9 


<0.5 


22 


14.4(0.6) 


5.4(0.8) 


2.0(0.4) 


3.6(0.5) 


7.1(0.5) 


0.6(0.2) 


<0.8 


<0.9 


1.3(0.3) 


<0.5 


23 


<1.9 


<2.2 


<1.3 


<1.6 


<1.5 


<0.5 


<0.8 


<0.9 


<0.9 


<0.5 


24 


<1.9 


<2.2 


<1.3 


<1.6 


1.8(0.5) 


<0.5 


<0.8 


<0.9 


<0.8 


<0.5 


25 


4.4(0.6) 


<2.2 


<1.3 


<1.6 


<1.5 


<0.5 


<0.8 


<0.9 


1.1(0.3) 


<0.5 


26 


<1.8 


<2.2 


<1.3 


<1.6 


2.9(0.5) 


<0.5 


<0.8 


<0.9 


<0.8 


<0.5 


27 


<1.8 


<2.2 


<1.3 


<1.5 


1.7(0.5) 


<0.5 


<0.7 


<0.9 


<0.8 


<0.5 


28 


<1.8 


<2.1 


<1.2 


<1.5 


<1.4 


<0.6 


<0.7 


<0.9 


<0.8 


<0.5 


29 


<1.9 


<2.2 


<1.4 


<1.5 


<1.6 


<0.9 


<0.9 


<1.0 


1.2(0.3) 


<0.8 


30 


<1.9 


<2.2 


<1.4 


<1.5 


<1.6 


<0.9 


<0.9 


<1.0 


<1.0 


<0.8 


31 


9.5(0.6) 


2.9(0.7) 


2.2(0.4) 


6.4(0.5) 


7.1(0.5) 


1.6(0.2) 


<0.8 


<0.9 


<0.8 


<0.6 


32 


2.7(0.6) 


<2.1 


<1.2 


2.1(0.5) 


<1.4 


<0.7 


<0.8 


<0.9 


<0.8 


<0.6 


33 


43.2(0.7) 


22.7(0.8) 


11.4(0.5) 


25.6(0.6) 


31.8(0.5) 


5.3(0.2) 


1.2(0.3) 


3.3(0.3) 


1.3(0.3) 


<0.5 


Line Width' 


; 235(3) 


202(8) 


277(9) 


132(4) 


253(4) 


546(22) 


357(37) 


305(27) 


361(30) 


551(75) 



1 Intensities given are the peak in mK. The format is intensity with the 1 sigma uncertainty in paratheses. For for non-detections, the 3 sigma 
upper limit on the line intensity is given. 

2 Line widths arc in km s — 1 . The format is line width with the 1 sigma uncertainty in paratheses. 



TABLE 4 

Assumed Molecular Abundance Relative to H2 



Molecular Species Ratio - GMCs 1 Ratio - NGC 253 2 



13 co 


2xl0~ 6 




HCN 


9xl0~ 9 


5.0xl0" 9 


HCO+ 


2xl0~ 9 


1.6xl0- 9 


HNC 


2xl0~ 9 


l.OxlO" 9 


cs 


6xl0~ 9 


6.2XHT 9 


N 2 H+ 


3xl0~ 10 




HC 3 N 


4xl0~ 10 


6.3xl0~ 9 


SO 


lxl0~ 9 


1.3xl0~ 9 


CH30H 


8xl0~ 9 


1.3xl0~ 9 


C 3 H 2 


5x10" 10 


5.0X10" 10 


HNCO 


lxl0~ 9 


1.6xl0" 9 



1 Abundances from lBergin et ahl (1997) 
Abundances from lMartin et al.l ^OOg ) 
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TABLE 5 
Mass and Beam Filling Fraction 



Galaxy 


D s 


^area 


Mass from 13 CO 


Sco 1 


Mass from CO 


Dense Gas Fraction 




kpc 




M Q 


Jy km s 1 


M Q 




NGC 253 


0.9 


0.26 


2.7xl0 8 


l.lxlO 4 


1.2xl0 9 


0.23 


Maffei 2 


0.6 


0.11 


6.2xl0 7 








NGC 1068 


3.3 


0.08 


1.1x10 s 


2.3x10 s 


1.8xl0 9 


0.61 


IC 342 


0.8 


0.08 


6.2xl0 7 


1.5xl0 3 


1.2xl0 8 


0.52 


M82 


0.8 


0.19 


1.7xl0 8 


6.7xl0 3 


5.8x10 s 


0.29 


NGC 3079 


4.5 


0.07 


2.0xl0 9 


1.6xl0 3 


4.4xl0 9 


0.45 


NGC 3690 


10.4 


0.01 


1.5xl0 9 


4.5xl0 2 


6.5xl0 9 


0.23 


NGC 4258 


1.8 


0.02 


1.1x10 s 


3.0xl0 2 


1.3x10 s 


0.85 


Arp 220 


17.5 


0.01 


4.7xl0 9 


4.0xl0 2 


1.6xl0 10 


0.29 


NGC 6240 


23.2 


<0.01 


<4.9xl0 9 


2.4xl0 2 


1.7xl0 10 


<0.29 



1 From lYoung et al.l ( 199")) 



